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ABSTRACT: We extend previous work to consider the effect of the soliton on the 
density of states and conductivity of quasi-one-dimensional Peierls systems with quantum 
lattice fluctuations, modeled by a random static disorder. Two features have been verified 
over an order of magnitude variation in the disorder. (1) The soliton density of states 
and the leading edges of both the soliton-to-band and the band-to-band conductivities 
have universal scaling forms. (2) The soliton-to-band conductivity has the remarkable 
feature that the leading edge is accurately predicted by the joint density of states while the 
trailing edge tracks the rigid-lattice conductivity. Or, in other words, disorder dominates 
the leading edge, while matrix element effects are predominant for the trailing edge. 

PACS: 71.20.Hk, 71.45. Lr, 71.55.Jv, 73.20.Dx 



I. INTRODUCTION 



This paper is the next logical step in working out the consequences of a simple observa- 
tion for most quasi-one-dimensional systems: the lattice zero-point motion is comparable to 
the lattice dimerization associated with the Peierls ground state. 1 The treatment of quan- 
tum lattice fluctuations about the dimerized ground state has to date resisted theoretical 
solution. Below we describe what has been accomplished. 

The first step was to model the quantum fluctuations with a static random potential. 1 
This proved possible, for temperatures below the transition temperature, since (i) the 
characteristic phonon energies are less than the Peierls energy gap parameter A and (ii) the 
phonons are then effectively dispersionless. At the mathematical level, the gap parameter 
is supplemented by a complex random variable SA(x) with a zero average and 

(6A(x)6A*(y)) _ 

^2 = V o{(x-y)/€o)- (1-1) 

By measuring length in terms of the coherence length £ G = hvp/A, we can characterize 
the static disorder by the dimensionless number rj. 

This model can be fitted into the standard continuum form for the electronic part 
of the Peierls- Frolich Hamiltonian 2 ' 3 (see Section III). The density of states p 'disorder f° r 
the resulting problem can be solved analytically, either for the case of an incommensurate 
lattice distortion 4 or for the commensurate case of a half-filled band. 5 Both cases give 
similar results for weak disorder but differ qualitatively for the experimentally uninteresting 
case of strong disorder (77 > 1). 

In comparison the density of states for the rigid-lattice 

1 1 \E\ 

Prigid-lattice( E ) = ^jT y /^~T ^~K'^ E ^ ~~ 

with a gap for energies \E\ < A has two important features: (1) the square-root singularity 
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at E = ±A is smeared out and (2) the density of states is now finite for all energies and 
has a substantial subgap tail. 

The second step was to extend model (described in detail in Section III) calculations to 
compute the optical conductivity. For the rigid-lattice the gap in the density of states leads 
naturally to a gap in the optical conductivity until Tiu = 2A where there is a square-root 
singularity. In contrast all experiments show a strongly rounded peak about 2 A, and it 
was known 6 that this shape could not be fitted by just rounding the rigid-lattice results. 

A complex analytical and numerical calculation 7 leads to a rounded peak conductivity 
which resembled existing experimental data. A striking feature was that the conductivity 
below the absorption peak could be fitted with a universal form 8 which could be used to 
deduce the Fermi velocity and r\ from the data. This work has encouraged more data from 
the experimentalists now able to fit their data.^' 1 ^' 11 

The third step - in this paper - is to compute the effect of solitons on the optical 
conductivity of quasi-one-dimensional compounds. Since solitons are the most studied and 
discussed feature of these materials, this work is long overdue. As Section III shows a 
considerable enhancement of the techniques developed before 7 was required. The next 
section summarizes the principal conclusions. More extensive results are given in Section 
IV followed by conclusions. 

II. SUMMARY OF RESULTS 

A. Review of soliton absorption in rigid-lattice system 

We use this section both to define notation and to describe the physics of the rigid 
lattice in a form appropriate for comparing with our results for disordered systems. As 
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Heeger et al. 12 observed a midgap absorption is a universal feature of actual Peierls systems, 
since solitons can be created by the formation of the material (e.g., in isomerization of cis 
to trans polyacetylene) or, more usually, by doping the material. 

The soliton is located midgap with a wavefunction proportional to 1/ cosh(a;/£ ) cen- 
tered in a CDW strand (— L < x < L) [£ = hvp/A]. The local density of states is 

■soliton 



sonion (j? \ _ -v"/ 

Prigid-lattice^iX) ^ 2 ^2 {x / & 



^ ) 

+ 4 V^ - A 2 ~ M2cosh^/u) ^' " A) " (2 ' 1} 
The density of states in the absence of the soliton is given by omitting both terms with 

cosh 2 (;r/£ ) in (2.1). (Note: f^dx sech 2 (a;)/2 = 1.) 

Using this form we could compute the local joint density of states which would indicate 
an absorption due to the soliton starting at A and localized in the vicinity of the soliton. 
More accurately one can compute the conductivity from the golden rule. 13 In particular 
the matrix element is proportional to the overlap of the continuum state and the soliton 
wavefunction: 

f ^ e tkx sech(x/Q = Trsech (ttA;£ /2) = V2^ k . (2.2) 
J £o 

Using the excitation spectrum above the gap E(k) = a/ A 2 + (%vp) 2 k 2 we write 

oc J dk\i) k \ 2 5{Tiu - E k ) (2.3) 

(2.4) 



soliton ( \ 
G rigid-latticed > 



soliton— to— band 

ioe 2 toK l 1 



2L h A 2 ^(%uj) 2 /^ 2 - 1 cosh 2 (7r V(^) 2 /A 2 - 1/2) ' 
In the last line we give a dimensionally correct form by including the "cross-sectional 

area" A of the one-dimensional strand. For reference we note that the integrated weight 

of the soliton conductivity is 2.83 (e 2 /h)(^ 2 / (2LA)). Most importantly, the conductivity 

is proportional to the density of solitons £ G /(2L). 

The dashed line in Figure 1 shows the soliton-to-band conductivity. We note for 
reference there is no analytic form for the band-to-band conductivity. Following many 
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workers 14 we use an approximate form for this conductivity 

n(l]Y oUton = e^_o (2A\ 2 2A / _ & 

rigid-lattice band _ to _ band h A \hu J ^/(htv)' 2 - (2A) 2 V V 

Equations (2.4) and (2.5) will be used for comparion with our results. 



B. Sample result for absorption in disordered system 



As discussed in the introduction, the effect of the zero-point motion is approximated 
by a static disorder which is characterized by a dimensionless parameter rj. In section III 
we will give the details of the procedure by which we compute the disordered case. The 
important point is that the calculation gives results which could not be anticipated by just 
smearing the rigid-lattice case - that is, with no disorder. As noted in the introduction, 
in the calculation of CDW conductivity, 7 such a convolution gives a large low frequency 
conductivity that is not observed and does not even reproduce the qualitative shape of the 
frequency-dependent conductivity seen experimentally. 

Figure 2a shows the soliton density of states for rigid-lattice and the disordered-lattice 
with a range of 77's. In Section IV we will discuss curves for a variety of 77 and show 
how the soliton peak can be scaled to deduce the dependence of the width on rj. Here we 
concentrate on the main features. The most important feature is the delta function midgap 
state (2.1) is broadened. This broadening has two consequences: (1) the soliton-to-band 
absorption will certainly be broadened on the leading edge and (2) the conductivity will 
be affected by the doping or Fermi level. When the Fermi level is tied to the middle of 
the soliton peak, only filled soliton states below the Fermi level can have transitions to 
the empty conduction states. On the other hand, if the Fermi level lies wholly above the 
soliton band, then the whole soliton band can have transitions to the conduction band, 
thus leading to a broader absorption than when the Fermi level is pinned to the midgap. 
This effect is real and should be observable. It is discussed in detail in Section IV. B. 
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In Figure 1 the solid line shows the conductivity for the disordered lattice with the 
same value of r\ = 0.01 as for the narrowest density of states in Fig. 2. Comparing with the 
dashed curve for the rigid-lattice we see that the leading edge is broadened as expected. 
In contrast the trailing edge is not appreciably broadened relative to the rigid-lattice case. 
We will return to this in the next subsection. In section IV there are more detailed 
comparisons. 

C. Comparison of rigid-lattice and disordered case 

We consider one more curve in Figure 1, the dotted line for the the joint density of 
states 

=JdE P S £ r °2er(E) P&SeriE + M- (2-6) 

Figure 3 shows the limits of the integration for transitions involving the soliton are 
such to permit two sets of transitions: (a) from the filled part of soliton band to empty 
conduction band and (b) from filled valence band to the empty part of soliton band. In 
the case where the Fermi level lies in the soliton band both transitions are possible and 
lead to a broader joint density of states. In the case where the Fermi level is wholly above 
(below) the soliton band only case a (b) is permitted and leads to a narrower joint density 
of states. 

The dotted curve in Figure 1 is the joint density of states for r\ = 0.01 with the Fermi 
level pinned to the midgap. The striking feature - nearly invisible - is that the leading 
edges of the joint density of states are nearly coincident with that of the full calculation. 
This is a common feature for all r\ and all chemical potentials (c.f. Section IV). This 
suggests that the matrix element is essentially constant for the leading edge. On the other 
hand, for the trailing edge the dotted line lies clearly above the full calculations. Indeed 
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as noted in the previous section the trailing edge tracks the rigid-lattice case which is 
dominated by the matrix element. In summary, 

1. The leading edge is dominated by the joint density of states for the disordered lattice. 

2. The trailing edge is dominated by full conductivity for the rigid lattice case. 

We have no simple explanation for relatively sudden switch on at fko = A of the matrix 
element effects and relatively sudden switch off of the disorder effects. But we can affirm 
that it has been seen in extensive numerical calculations over the parameter ranges 0.01 < 
rj < 0.5 and < E F < 0.5 A. 



We consider the standard continuum model of noninteracting electrons in one dimen- 
sion with a Peierls ground state. Our system is described by a Dirac-type equation for 
the wave functions ifti(x) and tp2( x ) f° r electrons moving with the Fermi velocity vp to 
the right and to the left, respectively, in the interval — L < x < L with a complex random 
potential 5A(x): 



In terms of the self-consistent disorder-averaged Peierls gap A in the ground state and 
the coherence length £ = hvp/A, the soliton configuration for the spatially- varying gap 
function A(x) is 



III. METHOD 




(3.1) 



(SA(x)) = (8A(x)8A(y)) = 0, (8A(x)5A*(y)) = r,A 2 8((x - y)/£ ). 



(3.2) 




(3.3) 
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We use G + (G ), the retarded (advanced) 2x2 matrix Green function, to compute 
all properties of interest. The density of states is 



p(E) = — f dx (--Im [Tr G + (x,x\E)}). 

2L J_ L 7T 



(3.4) 



The real part of the frequency-dependent conductivity for a strand of cross-sectional area 
A at low temperatures (ku> 3> kgT) is 



e 2 2 



E F 



dE — 



Ah nhu J-huj+Ep ' 2L J_ L 



dy 



L 



dx Re [j + (E + hu, E\x, y) - j ++ (E + fiw, E\x, y)] 



;++/ 



(3.5) 



The correlation function in (3.5) involves a trace of Green functions and the Pauli matrix 

°'3 



j ±± (E',E\x,y) = (hv F f(Tt [a 3 G ± (y,x\E')a 3 G ± (x,y\E)}) 



(3.6) 



Finally we compute the conductivity using the approximation (GG) (G)(G). 

Two linearly independent wave functions and V> satisfy the boundary conditions that 
ip(—L) and ip(L) do not diverge as L becomes large, whereas ip(L) and ip(—L) diverge. 
These i(j and -ip will be statistically independent in the limit L — > oo. 15 For x > y, 



G(x,y\E) 



4>l(x)ip 2 (y) ^i{x)ipi{y) 



(3.7) 



hv F (il>iil>2 - ihipi) \tp2(x)ip2(y) tp2(x)ipi(y) . 

where (ipifa ~ V^l) is a constant independent of x. We obtain G + (G~) by solving for 
wave functions with ImE a small positive (negative) number. Then the denominator can 
be Taylor-expanded and 



G + {x,y\E) 



oo 



V W pn(x,y)z n (x) p n {x,y)z n {x) 
hvF \^n(x,y)z n+ i(x) D n (x,y)z n+ i(x) 



(3.8) 



where 



C n (x,y) = y n (x) ^ 2 [ y \ , D n (x,y) = yn(x)^ 1 ' 



y n {x) 



ip 2 (x) 



z n (x) 
8 



^2 0*0 
ipl(x) 



(3.9) 



Next, we take the L — > oo limit and utilize the statistical independence of ip and ip 
to factor the average of the Green function into products of averages of C (D) and z. 
In general, the quantities C n (x,y) = (C n (x,y)) and D n (x,y) = (D n (x,y)) are functions 
of x and y, and = (yn(x)) and Zn(^) = (zn(x)) are functions of Thus the 

averaged Eq. (3.8) can be written in terms of z n (x), C n (x,y) and D n (x,y). The following 
equations for y n , z n and C n can be derived from Eq. (3.1) using standard Fokker-Planck 
methods, 16 subject to the conditions yo = zq = 1, dy n /dx\ x =±oc = dz n /dx\ x =±oo = 0, and 
C n (x, x) = y n { x )i which follow from Eq. (3.9) and the fact that the translational invariance 
has to be restored far from x = where the soliton is located: 



dy n E . A(x) . A(x) 2 

~dx U A Vn ~ m ~^~ Vn + 1 ~ m —^-Vn-l ~ m Vn, (3-10) 

t dz n . E . A(x) . A(x) 2 

~dx =2m A Zn ~ m ~^~ Z n+l ~ m —^~ z n-l ~ V n z n, (3.11) 

-??(n 2 + n + J j C„. (3.12) 



2 y 

The function D n satisfies the same equation as C n with the initial condition D n (x,x) = 

y n +i(x). 

We can relate y n and z n in a simple manner. In the soliton case with A(x) = —A(—x), 
we can easily prove that z n (x) = (—l) n y n (—x). In order to integrate Eq. (3.10) numerically, 
we choose a starting point (x = —L) sufficiently far from x = that the spatial variation 
of y n is negligible. Then y n satisfy a set of algebraic equations: 

2— + trin j y n - y n +\ - y n -\ = 0, (3.13) 

which we solve by truncating the number of equations. With the solutions of Eq. (3.13) 
as initial conditions, we integrate Eq. (3.10) from x = —L to +L. These results are used 
as initial conditions for Eq. (3.12). After we obtain y n (x), z n (x), C n (x,y) and D n (x,y), 
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and therefore G + (x,y\E) for x > y and E > 0, we get all other Green functions using the 
symmetry relations 





x\E) = 


G±(x,y\E) = [Gf 1 (x : 


y\E)T 




x\E) = 


G±(z,y\E) = [C% 1 (x, 


y\E)T 


G±(x, 


y\-E) 


= -G^ 2 (x,y\E), 




G±(x, 


v\-E) 


= G$ 1 (x,y\E). 





The key to our method is the fact that in disordered systems, y n , z n , C n , and D n all 
decay very rapidly as a function of n. To solve Eqs. (3.10-13) numerically, we truncate the 
number of equations by setting y n = z n = for n > N. The results presented here was 
obtained for N = 40, which was sufficiently large to ensure convergence. Our system size is 
20 coherence lengths, and therefore the soliton density is l/(20£ o ). In the numerical three- 
dimensional integration of Eq. (3.5), we used the mesh sizes of dE = 0.025 A, dx = 0.05£ o , 
and dy = 0.25£ o . Convergence of the data with respect to these parameters was also 
carefully checked. 



IV. RESULTS 

While the basic result was presented in Section II, here we present additional details 
on pertinent features. Among these are role of (a) disorder on the density of states, (b) 
chemical potential on the conductivity and (c) disorder on the soliton-induced conductivity. 

A. Disorder effect on density of states: scaling 



Figure 2a shows the density of states for values of the disorder from 0.01 < r\ < 0.08. 
The dotted line gives the density of states in the absence of disorder (Eq. (2.1)); of course 
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the delta function at E = cannot be exhibited. We see that the curves evolve smoothly 
with increasing 77. 

Figure 2b shows the scaling of the soliton band: p(E,rj)/p(E = 0, 77) versus E/T(rf) 
where the half-width r(?y) has to be deduced by fitting. Except for the largest r\ (0.15 and 
0.2) all curves lie on a universal shape which is nearly Gaussian, as evidenced by the fact 
that p(E = 0,?y)r(?7) varies less than one-half percent over the full 1] range. The insert in 
Figure 2b shows the variation of T(r]) with r\. The dotted line shows the efforts to do a 
similar scaling for the leading edge of the band gap. Within the scaling the two widths are 
very similar. 

B. Chemical potential: effect on conductivity 

Figure 4 shows how the chemical potential can broaden and shift the conductivity and, 
in the process, shows that the leading edge is predicted by the joint density of states while 
the trailing edge reflects the conductivity without disorder. There are two sets of three 
curves: for the chemical potential = and = 0.55 A. Each set has a solid line for the 
full conductivity, the dotted line for the joint density of states and dashed curve for the 
rigid-lattice conductivity. Figure 3 shows that the effect of moving the chemical potential 
from the middle of the soliton band to wholly above it would broaden the soliton-to-band 
conductivity. Measuring at half height the width from Tiu = A, we see the width is roughly 
doubled as the chemical potential moves from midgap to fully above the soliton band. As 
striking is how well the joint density of states mimics the calculated conductivity. This 
further demonstrates that the leading edge of the conductivity bears no effect of any energy 
dependent matrix element. Less impressive is the agreement of the trailing edge with the 
conductivity for the soliton without disorder (dashed line). 

Only space constraints prevent us from showing similar curves for a wide range of 
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rj. Of course the effect is less dramatic as the leading edge broadens but the agreement 
between pjDOS with the leading edge of the conductivity is common feature over the range 
0.01 < i] < 0.1. 



C. Disorder effect on conductivity 



Figure 5 illustrates the conductivity for four rj values (0.01, 0.02, 0.05 and 0.1). While 
the leading edges of the soliton-to-band and the band-to-band broaden as expected, the 
trailing edges of both are essentially unchanged, tracking the conductivity of the rigid- 
lattice result, shown by the dotted line. While the curve is for the chemical potential Ep = 
0.55 A, the results for the midgap chemical potential are similar. Figure 5b demonstrates 
that the leading edges of the soliton-to-band conductivities have a universal shape when 
rescaled. The insert shows the r\ dependence of the conductivity half- width T c (i]). A very 
similar scaling exists for the width of the leading edge of the band-to-band conductivity 
(dotted line in insert). 



V. DISCUSSION AND CONCLUSIONS 



First we should point out all the failings of the calculation. These are as follows: (1) 
We approximate (GG) = (G) (G). While in principle we know how to overcome this, the set 
of differential equations to be solved is sufficiently complex that we have not undertaken 
it. (2) We have used the rigid-lattice results for the variation in A(x) for the soliton. 
Once we have introduced disorder we should self-consistently solve the soliton energy and 
wavefunction. This admirable task is also numerically possible, but we have not attempted 
it. (3) We approximate the quantum fluctuation of the lattice motion by a random static 
disorder. We believe this is a good approximation in the low-temperature regime but have 
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no hard estimate of the error. (4) We have ignored any electron-electron contribution to 
any of the density of states or conductivities. We would dearly love to consider a tractable 
model that includes them. 

There is another problem with our startling result that there is no energy- dependent 
matrix-element effect on the leading edge of the soliton-to-band conductivity. It is solely 
the results of numerical calculation. We cannot yet produce an analytic or physical argu- 
ment. Given that we have been working of the ramifications of the large-zero point motion 
since late 1991, we do not see a way to this clearly desirable deeper understanding. We 
welcome the advice and comments of colleagues in our continuing efforts. 

Nonetheless we think there is a very interesting result - the turn off of the disorder 
effects and and the turn on of the matrix-element effects at %uj ~ A - which is aptly 
summarized at the end of Section II. We would strongly encourage experimentalists to 
produce better data on soliton-to-band conductivity on better samples. We have been 
impressed by the high-quality data 9 ' 10 ' 11 our earlier work 7 on the optical conductivity of 
the Peierls state induced. 

Finally we suggest one way that offers an approximate approach to a parameter-free 
analysis. In our earlier calculation 7 for the conductivity of the Peierls state, it was possible 
to extract rj from the scaled width of the leading edge of the conductivity. In this paper 
we have extended the calculations to smaller ry's (0.01-0.1) than in our earlier work (0.1- 
1.0). The widths for the leading edge of the band-to-band conductivity with and without a 
soliton roughly agree. Accordingly the measured band-to-band width (r(?y)/A = 0.83?y 0-60 ) 
could be used to extract an r\ that in turn could predict the soliton-to-band conductivity. 

This work was supported by the KOSEF (Grant No. 951-0209-030-1), Ajou University 
Research Fund, and the DOE - Basic Energy Sciences, Division of Materials Sciences. 
Cray time on the Ohio Supercomputer Center is deeply appreciated. 
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FIGURE CAPTIONS 



Fig. 1. Frequency-dependent conductivity for system of one soliton in a Peierls state 
for a one-dimensional chain of length 20£ o . Solid line corresponds to disorder parameter 
rj = 0.01, while dashed line to rigid-lattice (77 = 0). Dotted curve is the joint density of 
states for 77 = 0.01, scaled to match conductivity peak at E = A. Note leading edge of 
conductivity is indistinguishable from joint density of states while trailing edge is similar 
to rigid-lattice case. These two features are seen for all disorder parameters and for all 
chemical potentials (doping levels). In this case, chemical potential lies in the center of 
soliton band. Small structure at fiu = is a numerical artifact. 

Fig. 2. Density of states for system of one soliton centered in 20£ o chain of Peierls state 
for disorder parameters 77 = (dotted line), 0.01, 0.02, 0.04 and 0.08. The widths of the 
soliton and band density of states increase with r\. In (b) the scaled soliton density of states 
p(E,rj)/p(E = 0,77) is plotted versus E/T(rj), thus demonstrating the universal shape for 
77 < 0.1. Only large disorder (77 = 0.15, 0.2) cases break away from universal shape for 
E/T(rj) > 2. The insert shows the 77 dependence of T(?7)/A for both soliton band (solid 
line: 0.465 ?7°- 49 ) and leading band edge (dotted line). 

Fig. 3. Schematic plots of the transitions involved in the joint density of states for the 
E F = and 0.55 A. The symbols s-b and b-b correspond to soliton-to-band and band-to- 
band transitions. Since the Ep = soliton-to-band transition covers only half the soliton 
band, that joint density of states will be roughly half as wide as the Ep = 0.55 A curve 
as reflected in Fig. 3. 

Fig. 4. Fermi-energy dependence of the soliton-to-band conductivity for disorder parameter 
77 = 0.01. For each solid line, corresponding to Ep = and 0.55 A, there is a nearly 
coincident joint density of states for the same chemical potential. This coincidence is seen 
for the full range of 77 values. The rigid-lattice conductivity (dashed line) does not track 
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the trailing edge of the disordered-lattice conductivity with the same fidelity as the joint 
density of states does the leading edge. 

Fig. 5. Conductivity for disorder parameters r\ = (dotted line), 0.01, 0.02, 0.05 and 0.1; 
all for Ep = 0.55 A. In (a) note that rigid-lattice case (dotted) closely tracks the trailing 
edges of all the disordered-lattice cases. In (b) the scaled soliton-to-band conductivity 
(j(u;)/(Tp ea fc(?7) is plotted versus (cu — ^ pea k)/^c(v)j thus demonstrating the universal shape 
for the leading edge. The insert shows the r\ dependence of r c (?y)/A for both soliton-to- 
band (solid line: 0.685 ^°- 513 ) and band-to-band (dotted line). 



17 




Fig. 1 




Fig. 2 




Fig. 4 



12 3 

oj/A 




-2-10 1 2 

[a)-U)peak(r])]/^c(r)) 



Fig. 5 



